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The bound kaon approach to the strangeness in the Skyrme model is applied to investigating the 
possibility of deeply bound ppK~ states. We describe the ppK~ system as two-Skyrmion around 
which a kaon field fluctuates. Each Skyrmion is rotated in the space of SU(2) collective coordinate. 
The rotational motions are quantized to be projected onto the spin-singlet proton-proton state. 
We derive the equation of motion for the kaon in the background field of two Skyrmions at fixed 
positions. From the numerical solution of the equation of motion, it is found that the energy of K~ 
can be considerably small, and that the distribution of K~~ shows molecular nature of the ppK" 
system. For this deep binding, the Wess-Zumino-Witten term plays an important role. The total 
energy of the ppK~ system is estimated in the Born-Oppenheimer approximation. The binding 
energy of the ppK~ state is B.E. ~ 126 MeV. The mean square radius of the pp subsystem is 
\f{rlp) ^ 1.6 fm. 
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I. INTRODUCTION 



For recent years, lots of theoretical or experimental efforts to explore the possibility of nuclear .ftf-bound states 
[l| have been made. Although no firm evidences to show their existence are known up to now, there is one result 
reported by FINUDA collaboration which may suggest the existence of the lightest nuclear J^-bound state, ppK~ . 
In the experiment at DA$NE, it was observed that A and p from K~ absorption on 6 Li, 7 Li and 12 C at rest has a 
strong back-to-back correlation, and that the invariant mass spectrum of A and p shows a peak. The collaboration 
advocates that the observation can be interpreted as a signal of the formation of deeply bound ppK~ state, whose 
binding energy is 115 MeV and the width is 67 MeV. 

This experiment is motivated by the idea proposed by Akaishi and Yamazaki (AY) [l[ suggesting the existence 
of deeply bound K nuclei. It is based on the assumption that A(1405) baryon is a KN bound state formed by the 
strong attraction in / = KN channel. One may then expect that when a K~ is injected in a nucleus it may attract 
surrounding p's to form a shrunk nucleus. The K~ is bound deeply so that an absorption reaction, K~p — > 7rE, is 
energetically closed, and accordingly, it can have a long life time in a nucleus. 

However, it has not yet been established that the peak observed by FINUDA really corresponds to the state proposed 
by AY. Magas et al. |3| claimed that the peak corresponds mostly to the process K~pp — * Ap followed by final-state 
interactions of the produced particles with the daughter nucleus. Even if we suppose the peak to be a ppK~ bound 
state, it is strange that it is much deeper than the original AY prediction [l|: the binding energy B.E. = 48 MeV and 
the width V = 61 MeV. Recently, two groups 0,13 performed Faddeev calculations. The authors in Ref. [i| obtained 
B.E. = 55- 70 MeV and T = 95 - 110 MeV. The result in Ref. are B.E. ~ 80MeV and T ~ 73MeV. They are at 
odds with both of the result by FINUDA and AY prediction. On the other hand, an attempt to describe the if-nuclei 
as A(1405)-heypernuclei has also been made in Ref.jfJ]. 

A new experiment are planned to be performed at J-PARC searching for the deeply bound ppK~ state by the 
missing-mass spectrum of the 3 He(in-flight K~, n) reaction, together with the invariant-mass spectra detecting all 
particles decaying from the ppK~ bound state. It would be naively expected that a clear signal for formation of 
kaonic nuclei appears in this measurement, since a lighter nucleus is chosen as a target. Indeed, it has been suggested 
theoretically that a distinct peak of the ppK~ bound state can be observed in the spectrum of the 3 He (in-flight K~ , 
n) reaction if some conditions for the ppK~ optical potential are satisfied Q • 
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Our interests in this paper are whether the deeply bound ppK~ state can be realized in the context of the topological 
soliton model of baryons, the Skyrme model [8|. For this purpose, we employ the bound kaon approach to the 
strangeness in the Skyrme model [3] , which describes hyperons as the bound states of an anti-kaon and a topological 
soliton of the pion field ("Skyrmion"). 

The Skyrme model is a low ener gy e ffective theory of QCD at large number of colors. In the limit of large number 
of colors, as was shown by t'Hooft [10|, QCD reduces to a theory of weakly interacting meson theory. The action of 
the Skyrme model is a chiral effective theory written in terms of the Nambu-Goldstone boson fields. Nucleons emerge 
as topological solitons of the SU(2)f sector in the Skyrme Lagrangian (Til. H3|. 

One of the way to introduce the strangeness to the model is to assume a kaon field fluctuating around the SU(2) 
Skyrmion (bound kaon approach). One finds the existence of bound states of the kaon and a Skyrmion, which can 
be identified to be hyperons. The lowest bound state has the quantum number I = 1, t = 1/2, where I is the orbital 
angular momentum of the kaon and t the combined angular momentum and isospin, T = L + I , respectively. The 
parity of the I = 1, t = 1/2 state is totally positive, which is assigned to positive parity hyperons. A notable feature 
is the presence of a bound state in the negative parity state, I = 0, t — 1/2, which lies above the I = 1, t = 1/2 
state. This state probably corresponds to A(1405) baryon. Whereas the constituent quark models have difficulties to 
describe A(1405), this approach predicts the static properties of A(1405) [l3[ as well as octet and decouplet baryons 
in good agreement with the empirical values. In addition, an interaction originating from the WZW term acts on the 
S = +1 state, e.g. pentaquark, repulsively, and the state is pushed away into the continuum, while the interaction 
acts attractively on the S = — 1 state, leading to the formation of the bound states. 

Thus, the bound kaon approach is a theory naturally describing both of the positive-parity hyperons and the lowest 
negative parity state, A(1405), on the same ground. It is also worth mentioning that this approach has no parameter 
once we adjust and e (for their definitions, see below) to fit the N and the A masses in SU(2)f sector. Therefore, 
the KN interaction, which is a key ingredient for the study of K nuclei, is unambiguously determined. In these 
respects, it is of great significance to investigate the issue of the exotic nuclei such as ppK~ , K~ppn, K~pnn and so 
on in the context of the bound kaon approach to the Skyrme model. 

We describe the ppK~ system as two-Skyrmion around which a kaon field fluctuates. Each Skyrmion is rotated in 
the space of collective coordinate and its rotational motion is quantized to be projected onto a relevant two-nucleon 
state. We adiabatically treat the nucleon-nucleon radial motion and derive the kaon's equation of motion when the 
position of the Skyrmions are fixed first. Then we obtain the energy of kaon as a function of the relative distance 
between the two Skyrmions, which tells us whether the kaon can be deeply bound to two-proton or not. For the 
existence of nuclear K bound states, it is necessary that K gains sufficiently large binding energy in nuclei. If such 
nuclei exist, the nuclear components rearrange themselves under the influence of the strong attraction in I = KN. 
Then the nuclear part in K nuclei is excited relative to the original nuclear system. Therefore, the energy gained by 
K must be large enough to compensate the energy loss of nuclear component and to deeply bind the total system. 
If a strong binding of K~ is possible, it is also an interesting subject how the mechanism responsible for the strong 
binding is explained in the solitonic picture of baryons. 

In our previous paper (l4| . we have presented the derivation of the kaon's equation of motion and its numerical 
solution. It was shown that a strong binding of K~ to pp can occur. Needless to say, it cannot be taken as an evidence 
for the actual deep binding of the ppK~ system until the two-proton dynamics is treated. In the present paper, we 
give a detailed description of our approach, and attempt to solve the dynamics of the radial motion of the protons 
under the strong attaractive force mediated by K~ , in addition to the ordinary nuclear force. Then we can estimate 
the binding energy of the total ppK~ system within the Born-Oppcnhcimer approximation. The possible structure 
of the ppK~ state is also discussed. 

The organization of this paper is as follows. In the second section, we derive the kaon's equation of motion, and 
show its numerical solutions in the third section. We solve the radial motion of the two-proton in the forth section. 
Discussion and summary are given in the fifth section. Full expression of the kaon Lagrangian and the usefull formulae 
for collective coordinate quantization are gathered in the appendices. 



II. DERIVATION OF THE KAON'S EQUATION OF MOTION 



Let us begin with showing how the K~ coupled to pp is described in the bound kaon approach to the Skyrme model. 
We consider two Skyrmions fixed at positions with the relative distance, R, and assume the presence of the kaon field 
fluctuating around the Skyrmions. The equation of motion for the kaon in the background field of the Skyrmions is 
then derived, from which we know the behavior of the K~ coupled to pp. 

The action of the Skyrme model is given by 

r= [ d A x(^tiid^d^u) + -^tr[d^uu^d„uu^ 2 \ + r SB + r wzw (i) 
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where U is the chiral SU(3) field built out of the eight Nambu-Goldstone bosons. Tsb is the symmetry breaking term 
fl3l ] given by 



SB 



Flm\ + 2F'im 



K""K 



4<S 



tr [U + U f - 2] 



24 



tr 



V2X 8 (U + C/t) 



Ft - Ft 



4cS 



-tr 



(1 - ^\ % ){Ud^d^U + U f df,Ud^U f ) 



(2) 



where JTVfK) an d Fjt(k) are the mass and the decay constant of the pion (kaon), respectively. The last term in Eq.(J2]) 
has a role to renormalize the kinetic energy term for the kaon, while the first two terms renormalize the mass term. 
Twzw is the Wess-Zumino-Witten anomaly action (l5j : 



wzw 



' 240tt 2 



< f xe ' xvot ^tT{U ] d lJ ,UU^d v UU^d ol UU^d UU^d 1 U) 



N c denotes the number of colors. 

We assume the following "product" ansatz for the chiral field representing KNN system, 

U = U(1)U K U(2), 



(3) 



(4) 



where U(l) and U{2) are the fields of the baryon number B = 1 SU(2) Skyrmions located at r(l) = r — R/2 and 
r(2) = r + R/2, respectively. Their explicit expressions are as follows, 



U(i) 



u(i) 




, u (i) = e iF ^ T - f '^, (» = 1,2), 



(5) 



where r(i) — \r(i)\, r(i) — r(i)/r(i), F(r) is the profile function of an isolated Skyrmion [12[ and r the Pauli matrices. 
Uk is the field carrying strangeness. Its form is 



Uk = exp 



. 2^2 / 
' F K 



K 



where K is the usual kaon isodoublet, 



K = 



K+ 
K° 



Each Skyrmion is rotated in the space of SU(2) collective coordinate, Ax or A 2 , as 

u(l) Aiu(l)A{, u{2) ^ A 2 u{2)A\. 



(6) 



(7) 



(8) 



By substituting the ansatz, Eq.((3]), with the replacement Eq.© into the action, Eq.((T]), we obtain the Lagrangian 
for the kaon field in the presence of the background B — 2 Skyrmion. After expanding the Lagrangian up to the 
seceond order in K, we obtain 

C = {D^ D^K - m 2 K K^ K - ^ K ^{d^Ul^UBB) + -i^tr^BB^L. d v U BB U BB } 2 ^ 
jL- \2{D lx KfD v KtT{A^A v ) + D^Ktt{d v U% B d v U BB ) - 6(D^[A" , A»]D V K 



iN 

-^B^D^K-iD^K], 

K 



where U B b represents the product of rotating solitons, 

U BB = A 1 u{l)A\A 2 u{2)A\. 
In Eq.©, the "covariant derivative" is defined by 



(9) 



(10) 



(11) 
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and 

V, = + i^(2)]/2, ^ = - i^(2)]/2, (12) 

where 

L„(l) = A 1 uHl)A\d fi [A 1 u(l)A\], i? M (2) = A 2 «(2)4^[A 2 «t(2)4]. (13) 

The last term in Eq. © comes from the WZW term and is the baryon number current given by 

Btl = e -^X^(ul B d,U BB ul B d a UBBU t BB d p U BB ). (14) 

Here we note that the Lagrangian, Eq.©, has the same form as that for B = 1 Skyrmion 0] except that the 
background Skyrmion is not the single B = 1 Skyrmion but the product of B = 1 Skyrmions, Eq. lfTU)) . It should be 
also noted that the KNN interaction is unambiguously determined as in Eq.© once the ansatz for U is given. 

We neglect the terms suppressed by 1/N C in Eq.©. Since the time-derivative of the collective coordinate is 
A lt2 ~ OiN- 1 ), we see A , V Q and Bj are 0(N^ r ) from their definitions, Eqs.lO and fli]). 

A ,V ~ OiN- 1 ), B^OiN- 1 ). (15) 

Then the Lagrangian for the kaon field up to 0(N®) terms reads as follows, 

2rt tX ( d 3 U BB d 3 U BB) 



C = (doK^doK 



2e 2 Fj. 



K^DjDjK - m\K^K 



\&k\ - tr^Ul^UBB) + -^tr^UsBUl^d^BBUls} 2 



2K^D j [D i Ktr{A j A i )} + ^^[^^(6,17^9,^)] - 6K*Dj([Ai, A 3 ]D % K) 



e 2F K{ - 3 L-.«-n"j"./j . 2 
i N 

-^[K^d K - (d K?K]B° + 0{N- 1 ). (16) 

K 

A comment is in order here. It should be noted that in the Lagrangian Eq. (|16[) Ubb must be the product not of 
the static solitons but of the rotating solitons, Eq. (TlT)]) . In the case of B = 1, the effect of rotation is suppressed in 
the limit of large N c since the collective coordinate enters the Lagrangian only through its time-derivative, which is 
0(N~ 1 ). As an example, lets us consider 

tr^f/^t/t], (17) 
which appears in the coefficient of K*K in Eq.©. For B = 1, under the replacement 

u -> AuAt (18) 

with A being a collective coordinate, Eq. (fT7|) reads 

tr^ud'V] -> tr [Ad lU A^ Ad 1 u^A^} + ©(A^ 1 ) = tr^ucM] + O^ 1 ). (19) 

Therefore, for B = 1, the soliton may be regarded as a static one at N c — > oo. On the other hand, in the present 
B = 2 case, since each Skyrmion is rotated independently, 



u -> Aitt(l)4^2tt(2)4- ( 20 ) 



Eq. p?)) reads as follows, 



tr[0 /1 u0 /i v t ] -» tr[9 4 (A 1 u(l)4^2w(2)A^)9 4 (A 2 w t (2)4^iw t (l)4)] 
= tr[(Ai5iu(l)4^2w(2)4 + ^(1)4^2^(2)4) 
x (^2^^(2)4^1^(1)4 + ^2^(2)4^1^^(1)4)1 

= tr[a,w(i)4^2w(2)a i w t (2)4^iw t (i) + 0^(2)0^(2) 

+fi| i u(l)e i « t (l) + w(l)4^29 i M(2)M t (2)4^ia l u t (l)] + 0{N~ X ). (21) 
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Thus the collective coordinates themselves, which are not suppressed by 1/N C , enters the Lagrangian. As long as using 
the ansatz, Eq.Q, the kaon inevitably couples not to static soliton but to rotating solitons in the limit of N c — ► oo. 
Accordingly, we consider the kaon under the background of the two-Skyrmion projected onto a spin-isospin state of 
two-nucleon. 

Our next task is to perform the collective coordinate quantization, and project the rotation of each Skyrmion onto 
the relevant spin-isospin state. This procedure is done as follows. First, we rewrite the Lagrangian, Eq. (|16| . in terms 
of the adjoint matrix defined by 

D ij (A)=tT(T i Ar j A^)/2 (22) 

with A being a collective coordinate. The result, which is quite lengthy, is displayed in Appendix [A"l The matrix 
Dij(A) is known to be represented by the rotation matrix of rank-1. On the other hand, the wave function of the 
nucleon in the space of collective coordinate is also expressed by rotation matrix [12j : 

(A\Nis,J 3 ) = ^(-l) /3+1/2 £ V iU(tt), (23) 

where 1% and J3 denote the third component of the isospin and that of the spin, respectively, and f2 the Euler angles. 
Then, the projection of the Skyrmions onto physical two nucleon states is performed by sandwiching the Lagrangian, 
Eq. p^)) . with two nucleon states and integrating the Euler angles, 

dfii J dn 2 {N{l)N{2)\C\N(l)N(2)), (24) 

where N(i) denotes the i-th nucleon. We assume that the proton-proton in the ppK~~ system is in spin-singlet, and 
project the rotational motion of the Skyrmions onto the spin-singlet proton-proton state. Thus we consider 



£ppK — 



J dSli J dn 2 {(pp) s=0 \C\{p P ) s=0 ), (25) 



where \(pp) s =o) is the wave function corresponding to the spin-singlet proton-proton state and is given by 

IW-o) - ^(|%i(l)iV W (2)) - 1^(1)1^^1(2))), (26) 



with Af being the normalization constant. Eq. (|2"5|) is the Lagrangian for the kaon coupled to two protons. Detailed 
description of the projection are shown in Appendix [Bl 

Now, we derive the equation of motion for the kaon from the Lagrangian, Eq. (|25[) . First, we average the direction 
of the line joining the two Skyrmions. To do that, we put R/2 = ((R/2) sin a cos f3, (R/2) sin a sin f3, {R/2) cos a) in 
the Lagrangian Eq. (f2l))) . and integrate the angles a and f3, 

£ PP K = — / da d/3sma£ ppK . (27) 
4?r Jo Jo 

This corresponds to assuming that the proton-proton system is in S'-wave. Then the background field becomes 
spherical, which allows us to set the kaon field as 

K(r,t)=k(r,t)Y lm {e,<f>) (28) 

with Yi m (9, 4>) the spherical harmonics. This ansatz, Eq. (|28| . is substituted into the Lagrangian, Eq. (|27|) . We perform 
the 9- and (^-integrations before taking the variation with respect to k(r, t). Up to this step, quite long and involved 
calculations are needed. Then the Euler-Lagrange equation for k(r,t) yields 

d 2 _ d 



f(r;R)^-2fX(r;R)j f -m 2 K -V eS (r;R,l)+0 fc(r,i) = 0, (29) 
where the operator O is defined as 



Cl (r;R)—+c 2 (r;R) — . (30) 
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In Eqs. (|29f and (|30|) . the coefficients, f(r;R) and X(r;R) are expressed as follows: 



f(r;R) = daj^ d^wna^l+^g-ixidjU^gdjUBB)^, (31) 



i r , .„ . -at, 



A(r;i?) = — / da / cZ/3 sina— ^ <B°) . (32) 



47T./o Jq F k 



Here (•••) means taking an expectation value with respect to the spin-singlet proton-proton state as in Eq. (|25l) . 
— V e s (r; R, I) and O correspond to the terms with and without spatial derivative in the following equation, respectively, 



y e ff(r; i?, 0+0 = / del d<f>smOYi m (d,(j>)- — da d/3sina 
■/o Jo 4?r 

^^ [Atr^ A,)] + ^[Dj-tr^t^Bft^Bfl)] - 6^([A, A,]A)}\ ^(9, $(33) 



ci and C2 in 0, and V c g(r; i?, Z) are read off from this equation. Their explicit expressions are quite lengthy and are 
not particularly instructive. Therefore we do not display them here. 
Let us expand the field k(r, t) in terms of its eigenmodes: 

fc(r,t) = [ fc «( r ) e ""'4 + fen(r-)e- i ^*6„] , (w n , £ n > 0), (34) 

where a„ and 6 n are the annihilation operators for the strangeness S = =Fl states, respectively. Substituting Eq. (|3"4"l) 
into Eq. (|2"9")l . we find the eigenmodes satisfy 

f(r-,R)Ljl-m 2 K -Vt\r-,u n ,R,l) + 6]k n (r) = 0, (35) 
/>; - m|- - l4 +) (r; u>„,R, I) + o] fc„(r) = 0. (36) 

Eqs. (l35|) and (f36|) are the equation of motions for S = — 1 and 5 = +1 states, respectively. In Eqs. (|35"|) and {36}, 
Vk- (r; cj, i?, Z) plays a role of potential term and can be separated into two terms, 

V K T) (r; u, R, I) = V^ T J W (r; w, i?) + V eS (r; R, I), (37) 

where V^^ w (r; w, R) originates from the WZW term and V e s(r; R, Z) from remaining terms in the Skyrme Lagrangian. 

Vwzw( r '> w ' R ) is s iven b y 

V™ w (r; uj, R) = T 2A(r; R)u. (38) 

Thus the WZW term acts on the negative (positive) strangeness states in attarctive (repulsive) way. k n (r) and k n (r) 
obey the following normalization conditions, 

Awr 2 dr [/(r; R)(u n + u n >) + 2A(r; R)] k n {r)k n ,{r) = <W, (39) 

Anr 2 dr [/>; R)(uj n + «„,) - 2A>; JZ)] fc„(r)fc n , (r) = <W- (40) 

III. NUMERICAL SOLUTION OF THE KAON EQUATION OF MOTION 

We solved numerically the equation of motion, Eq. (|35p . FigureQ] shows the obtained energy eigenvalue of K~, u>, 
as a function of the Skyrmion-Skyrmion relative distance, R. m^, and e were taken to be — 0, F n = 129 MeV 
and e = 5.45, which were adjusted to fit the masses of N and A [12|]. The kaon mass was taken to be rriK = 495 MeV. 
For the ratio, Fk/F^, we have examined two choices: Fk / 'F n = 1.00 [l9[ and the empirical value, Fk/F^ = 1.23. 
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FIG. 1: The energy eigenvalues of S- and P-wave K , ui, as functions of the proton-proton relative distance, R. Two choices 
of the ratio Fx/F n are examined: Fk/F^ = 1.00 and the experimental value, Fk / F n — 1.23. 
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TABLE I: Energy eigenvalue of S-wave K (uji=o) and that of P-wave (uii=i) for five cases of the proton-proton relative 
distance, R. Fk I 'F n is taken to be the empirical value, Fk/F^ — 1.23. 



Wc find that the lowest-lying mode is the S'-wave and that the P-wave lies above the S- wave. Note that this order is 
natural but different from the case of B — 1 , where the lowest-lying mode is P-wave, as mentioned in the introduction. 
We see that it is important to take into account the difference between Fk and F w . As was shown in Ref. [IH, by 
setting Fk/F v equal to the empirical vaue, Fk/F^ = 1.23, hyperon masse are well reproduced, while when we set 
Fk/F w = 1 they are overbound. The binding of K~ to pp is also weaker when taking the empirical value of Fk/F^. 

In Table [H the K~ energy eigenvalues for several values of R are displayed. Looking at the .S'-wave channel, the 
binding of the kaon is extremely strong for smaller distance, i.e. R < 1.0 fm. In this region, the repulsive nucleon- 
nucleon interaction dominates over the attractive interaction between K~ and pp. Nuclear matter in which average 
NN distance is R ~ 1.0 fm is expected to be realized in the core region of compact stars. Our result shows that 
K~ in such high density nuclear matter can be lighter than the pion. This might be somehow related with the kaon 
condensation [1 €31 ] which is considered to occur in compact stars. As R is increased, the binding becomes looser. 
However, at R — 2.0 fm, for instance, which is close to the average inter TV TV distance in normal nuclei, the binding 
is still deep: the binding energy is about 140 MeV (see Table [TJ> . 

Next, we consider the P-dependence of the i4T~-distribution. In Figure[2j we plot the distribution of K~ in S'-wave 
(kaon wave function k(r) normalized by the condition, Eq. (|3"fj)) ) for several cases of P, R = 1.5, 2.0, 2.5 and 3.0 fm. 
For comparison, the baryon number denisity given by 

(P°>= J dSli J dn 2 ((pp) s=0 \B°\(pp) s=Q ), (41) 

is also plotted. At P = 1.5 fm, it can be seen that K~ is localized in the narrow region between the two protons. 
At relatively larger separation, P ^ 2.0 fm, K~ has large probability to stay near the proton's respective positions, 
which is characteristic to molecular orbital states Il4j], 
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FIG. 2: S-wave K~ distribution (normalized wave function k(r)) and baryon number density (Eg. (|4ip ) for the relative distance 
of the two Syrmions, R = 1.5, 2.0, 2.5 and 3.0 fm. Both are in unit of eF n . The horizontal axis is the distance from the origin. 
Fk/F-k = 1.23 was chosen. The baryon number density is multiplied by a factor 10. 



IV. PROTON-PROTON RADIAL MOTION 



Now, we solve the dynamical problem of proton-proton radial motion under the strong attaractive potential mediated 
by K~ , in addition to the ordinary nucleon-nucleon potential. We assume that the radial motion of the two-proton 
is governed by the following Hamiltonian, 

H = T NN (R) + V NN (R) + V K [R) , (42) 

where Tn n is the NN kinetic energy term, 

T ™ - "JOT ( r2 m) ■ < 43 » 

Here, the nucleon is regarded as not a soliton with finite size but a point like particle with the mass Mm = 939 MeV 
and its motion is assumed to be non relativistic. Vnn(R) is the state-independent part of the NN potential obtained 
from the product of B = 1 Skyrmion We have neglected the contribution from the state-dependent part since 
they give a minor contribution compared with the state-independent part. The last term is the energy generated by 
the bound kaon in S-w&ve, 

V K {R)=Lu l=0 (R)-m K , (44) 

where uji =0 (R) is the S-wave kaon's energy obtained in the previous section. Vmn(R) + Vk{R) may be regarded as 
effective pp potential in the ppK~ system. In the left panel of Figure [31 we show the behavior of Vnn(R), Vk(R) and 
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FIG. 3: Potential terms in Eq, (l42[) as functions of the proton-proton relative distance, R. The upper curve is the proton- proton 
potential in the absence of K~ , Vnn{R)- The lower one represents the energy of K~ , Vk(R)- The middle one corresponds to 
their sum, Vnn(R) + Vk(R), the effective proton-proton potential in the ppK~ system. 



Fk/Ftt (T nn ) (MeV) (V NN ) (MeV) (V K ) (MeV) total (MeV) ^(r 2 NN ) (fm) 
1.00 50.4 97.3 -380.5 -232.7 1.46 

1.23 42.0 74.5 -239.2 -125.5 1.63 



TABLE II: Calculated total energy of the ppK bound state relative to 2Mjv + m,K and its decomposition. The root mean 
square radius of NN subsystem, \/Jr%~^), is also shown. For the definition of each component, see the text. 



their sum. Vnn(R) has repulsive core at short distances and one-pion exchange tail at long distances. Medium range 
attraction cannot be produced within the product ansatz of B = 1 Skyrmions [l8[ . On the other hand, the attractive 
force generated by bound kaon, Vk(R), is so strong that it overcomes the strongly repulsive Vnn(R)- As a result, 
the effective NN potential in the ppK~ system, Vnn(R) + Vk(R), is strongly attractive in the medium range. 
The energy of the ppK~ state relative to 2-Mjv + ft^Ki E, is obtained by solving the Schrodinger equation, 

H^ N =E^ N . (45) 

In Table HI1 we show the energy oippK~ bound state relative to 2Mjv+TOk and its decomposition into the NN kinetic 
energy, (Tnn), and the NN potential energy, (Vnn), and the kaon's energy, (Vk), obtained by solving Eq. (|45|) . The 
expectation value, (O), is defined by 

/>00 f-OQ 

(0) = J 4TrR 2 dR^N(R)*O^N{R)/ J 4:TTR 2 dR^N{R)*^N(R)- (46) 

The root mean square radius for NN subsystem, y/ (r NN ), is also shown. 

When the Fx/F^ is taken to be the experimental value, Ffc/F^ = 1.23, the binding energy of the ppK~ bound 
state is estimated to be 

B ppK - ~ 126 MeV. (47) 

Another fact worth noting is the smallness of the NN kinetic energy, (Tnn) ~ 40 MeV. As long as looking at this 
fact, the Born-Oppenheimer approximation seems to be not so poor. The root mean square radius of the two-nucleon 
subsystem is 

^1.6 fin, (48) 
which is significantly smaller than the average inter NN distance in normal nuclei. 
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(Fk/F^I.23) 

S-wave Q 

S-wave (WZW off) 

P-wave — t— 
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FIG. 4: The energy eigenvalue of K~ as a function of the proton-proton relative distance, R, compared with the case when 
the Wess-Zumino-Witten term is switched off. The ratio of the decay constant is taken to be Fk/F^ = 1.23. The bottom and 
upper curves are the energies of K~ in S- and P- wave, respectively. The middle one represents the energy of S-wave K~ when 
the Wess-Zumino-Witten term is switched off. P-wave K~ is unbound when the Wess-Zumino-Witten term is switched off. 



V. DISCUSSION AND SUMMARY 



Lets us consider the mechanism responsible for the strong binding of K~ to pp. This is turned out to be attributed 
to the WZW term. It is known that in B = 1 sector the existence of the WZW term leads to various important 
results. The WZW term is included in the action from the requirement that the effective theory written in terms of the 
Nambu-Goldstone bosons should reproduce the anomaly the fundamental theory possesses. The Skyrme Lagrangian 
without the term has a fictitious symmetry: an invariance under U «-> . This symmetry forbids processes changing 
even-oddness of meson number, e.g. K + K~ — > 7r + 7r~7r°. The WZW term breaks this extra symmetry [l5j . In 
addition, for odd N c , the WZW term ensures that the quantized Skyrmion has a half-odd spin and thus behaves 
as a fermion. The effect of the WZW term goes beyond these rules. In the bound kaon approach, the interaction 
Lagrangian of the kaon and the nucleon which comes from the WZW term has the form of 

7 N 

Avzw = --^[K^d K - (d K)^K}B a . (49) 

This is nothing but the so called Tomozawa- Weinberg term. This term gives an effective attractive contribution to 
negative strangeness states which is crucial for obtaining the correct values of the masses of ground state hyperons. In 
particular, without the term, the S*-wave bound state of a kaon and a Skyrmion, which corresponds to A(1405), does 
not exist. Also in the present B = 2 case, the role of the WZW term is revealed to be important. In the equation of 
motion, Eq. ([35|) . there exist two terms, Vwzw ano - ^c ff m Ea. ([37|) . which effectively play a role of the potential acting 
on the kaon. Among them, it is Vwzw ^ na * originates from the WZW term. In order to see the effects of the WZW 
term, we switched off V^zw ana - calculated the kaon's energy. The result for the S-wave K~ is shown in Figure 0] 
One observes that the WZW term additionally gives a substantial attractive contribution to the binding of the kaon. 
Here, we should note that V^zw ^ s stronger than that for B = 1 since the interaction Eq. (f4"9")) is proportional to the 
baryon number density. It is therefore quite a natural result that K~ is bound to two-proton more deeply than to one 
proton. However, this does not necessarily mean that the larger the baryon number B becomes the deeper the binding 
of the kaon is, since V^zw does not necessarily become stronger with the increasing B: V^zw IS proportional to the 
kaon's energy (see Eq.flSS])). 

The X _ -distribution shown in Figured] suggests that the ppK~ state is a molecular orbital state. This observation 
is quite natural in the following sense. The two protons in the ppK~ system should keep some distance so as to 
avoid the repulsive core of the nuclear potential. On the other hand, for the two-proton with finite separation, it can 
be shown that the potential acting on the kaon, vjf \r: to, R, I) in Eq. ([57|) . is a double-well potential which is most 
attractive at the proton's respective position. It is the molecular orbit that the kaon occupies under such a double- 
well potential. If the ppK~ is really a molecular orbital state, it is plausible that the binding of K~ to two-proton is 
stronger than to one proton since K~ experiences the strong attraction from the two protons without increase of the 
kinetic energy. 
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Finally, an important comment is in order. As was discussed in the section [TT1 the kaon is inevitably bound to 
rotating solitons in the limit of N c — + oo, as far as we employ the product ansatz, Eq.(H]), and rotate the solitons 
independently. Therefore we have projected the Skyrmion rotation onto a relevant spin-isospin state of the two- 
nucleon, spin-singlet proton-proton state. Thus the system which we have considered is purely a ppK~ bound state. 
On the other hand, in the original bound kaon approach, hyperons do not necessarily correspond to KN bound states 
but quantized states of the bound system of a SU(2)f soliton and kaon fluctuating around the soliton. It would be 
appropriate to investigate the ppK~ states by extending such a picture of the hyperon to B = 2 and S = — 1 systems, 
which is left for our future work. 

In summary, we have applied the Skyrme model to a study of the lightest if- nuclear bound state, ppK~ . We have 
derived the equation of motion for the kaon coupled to two-proton at fixed position. The two-proton is expressed 
by two-Skyrmion whose rotational motion in the space of collective coordinate is quantized and projected onto spin- 
singlet proton-proton state. Numerical solution of the equation of motion shows that K~ can be strongly bound even 
for relatively large inter proton-proton distances and that K~ is in a molecular orbital state. Next, we have solved the 
two-proton radial motion by assuming that the protons move under the strong attractive potential generated by K~ in 
addition to the ordinary NN potential. Then we have found that ppK~ state can be realized as a very deeply bound 
and compact state, whose binding energy is B ppK ~ ~ 126 MeV and the mean inter NN distance is yj {r 2 NN ) — 1.6 
fm. The obtained value of the binding energy is surprisingly close to the experimental result obtained by FINUDA 

^exp. 

ppK~ 

this agreement. 



collaboraton, B cxp K _ =115 MeV. However, considering the crudeness of our treatment, we are not allowed to satisfy 
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APPENDIX A: FULL EXPRESSION OF THE KAON LAGRANGIAN 

In this appendix, we rewrite the Lgrangian, Eq. (|16[) . in terms of the adjoint matrix, Eq. (|22p . The adjoint marix is 
included in Ubb, Vj and Aj. Vj and Aj can be rewritten using the adjoint matrix as follows, 

Vj = iTm^D^A^L^l) + D mi (A 2 )R)(2)], (Al) 
Aj = iT m i[£> mi (Ai)Lj(l)-D mi (A a )J2j.(2)]. (A2) 



Here L* (l) and i?*-(2) are defined as 

4« 



j 

1 „ 1 



S(l) jiV ^-C(l)S(l) + f (l) 3 f (l)^'(l) + ejik f(l) k — 5(1) 2 , (A3) 
-Lc(2) S(2) _ f(2)/ (2)^(2) + e iifc f(2) fe -l . 



R){2) = -5{2) ] —C{2)S{2) - f(2)/(2) 4 F'(2) + e yik r{2) k — S{2) 2 , (A4) 



where 



6(a)ji = Sji - f(a)jf(a)i, (A5) 
F(a)=F(r(a)), F '(a) = ^^, (A6) 

C{a) = cosF(a), S(a) = sinF(a), (a = 1,2). (A7) 

Now we use these equations to rewrite the Lagrangian in terms of the adjoint matrix and Skyrmion profile function. 
Let us show the result for each term in Eq. p^)) . 

1. (9 ^) t 9o^-term 

ti(djUg B djUBB) in the first term of Eq. (TT6|) is expressed as follows, 

tiid 3 U BB d,U BB ) = 2 [4(1)4(1) + 4(2)^(2) - D ik (4A 2 )2LUl)R*(2)] . (A8) 
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2. K^K-term 



The first term in the curly bracket of the fourth term is given in Eq. (|A8p . The second term in the curly bracket is 
rewritten as follows, 

-Jpz tr [dp U B b U bb , d v U B b U bb ] 2 

= ^{[(^(^^(l) + ^(2)^(2))(4(1)^(1) + 4(2)^(2)) 
-(4(l)i}(l) + ^(2)i2}(2)) 2 ] 

+D fc „(4 J 4 2 )[-2(^(l)^(2) + L k i {l)BJy{2)){L\{l)L\{l) + R\(2)R l 3 (2)) 
+4^(1)^(2)^(1)^(1) + R l m (2)R l m {2))] 

+DUAM)D lm {A\AML k j {m-m + L k {l)R^2)){L\{l)Rf{2) + L]{\)R™{2)) 
-4^(1)^(2)^(1)^(2)]}. (A9) 



3. Terms containing DjK 

In this subsection, we display the results for the terms containing spatial "covariant derivative", DjK. The second 
term in Eq. (|16|) is rewritten as 

tfDjDiK = \{L){\)L){1) + R){2)R){2))] + iT m D mi {M)\[d 3 L){\) + L){l)d 3 ] 

+iT m D mi {A2)\[d R){2) + R)(2)d 3 ] + D lk {A\A 2 )[^L){l)R k J (2)]}K. (AlO) 

The three terms in the curly bracket of the fifth term have very lengthy expressions. 
The first term in the bracket: 

2K^D j [D l Kti{A J A i )] 

= 2K^[-^d j (L\(l)Ll(l) + R l j (2)R l i (2))d i 

+~(L?(1)L?(1) + i$(2)J2? (2))(4(1)Z&1) + R){2)R\{2))] 

+D lk {A\A 2 ) l -d 3 (L l 3 (\)R>y{2) + L\{l)R k j {2))d l 

-D^A\A 2 )D lk {A\A 2 )\{LV{l)R?{2) + L?(l)J^(2))(4(l)iZ*(2) + L\(l)R k (2)) 
-^r„^ ^^ (A 1 )i[a J Lf(l)(4(l)L ^ ^ (l) + ^(2)^(2)) + (Lj(l)L{(l) + ^(2)^(2))i?(l)9 i ] 
-ir^^aji^JJ? (2)(4(l)Li(l) + i2j(2)i&2)) + (L$(l)^(l) + ^(2)^(2)^(2)$] 

+ir„D n8 (^)A fc (4^)|[^iI(l)(4(l)^(2) + ^(l)^(2)) 

+(4(1)^(2) + Li(l)4(2))L?(l)5 i ] 
+ir n I? ri9 (A 2 )A fe (4^)|[5^!(2)(4(l)^(2) + I4(1)J#(2)) 

+ (4(l)i?f(2) + Li(l)4(2))i?|(2)^}K (All) 

The second term: 

iift^.^.^tr^f/^^t/Bs)] 

= Xt{[^(L« (l)L^(l) +J R« (2)ift(2))fi, - \{L){\)L){1) + R)(2)R)(2)) 2 } 

-D nk (A\A i )2d j L^{l)B^(2)d j 
+D <J (4A 2 )D nfc (AlA 2 )i5(l)i^.(2)i»(lX(2) 

+i7iAi(A 1 )i[L5(l)(L» (l)L^(l) + iC(2)iC(2))d, 
+5,(^(1)^(1) + i£(2)i£(2))L*(l)] 
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+inD u {A 2 )-[R)(2){L n m (l)L n m (l) + i£(2)i£(2))0, 

+9,(^(1)^(1) + JC(2)iC(2))i?}(2)] 
+ir i Ai(A)^(4^)2[-4(l)^(lX(2)^ - ^"(1)4^(2)4(1)] 

+*r i A 4 (A 2 )A lfe (A t 1 A 2 )2h^(2)L^(l)i? ? fe „(2)^ - 9,^(1)^(2)^(2)]}^. (A12) 

The third term: 

= -3Jft{-i[(4(l)Lj(l) - i?$(2)^(2))(L*(l)L*(l) - J2*(2)J2*(2)) 
-(4(1)4(1) - 4(2)i?f (2))(L< (1)L<(1) - 4(2)4(2))] 

-~[4(1)(4(1)4(1) - 4(2)4(2)) - 4(1)(4(1)4(1) - 4(2)4(2))]^ 

+^•[4(1X4(1)4(1) - R\{2)R]{2)) - L){1){L\{\)L\{1) - 4(2)4(2))]] 
~iT m D mk (A 2 ) [e knl d, Rf (2)4 {2)8, 

+^(2x4.(1)4(1) - 4(2)4(2)) - 4(2)(4(i)4( 1 ) - R\m^)m 
-is i [^(2)(4(i)4(i) - 4(2)4( 2 )) - R-(m\m\(i) 4(2)4(2))]] 

+iT m D mk (A 1 )D lr (A\A 2 )[e kln d j L](l)Rl(2)d i 

-i4(i)(4(i)4(2) - 4(W (2))^ - ^4(i)(4(W(2) - 4(i)4( 2 ))] 

+ir ro £> mfc (A 2 )Ar(4^2)[efem5 , i 4(l)4(2)a i 

+^4(2)(4(i)4(2) - 4(1)^(2))^ + ia,4(2)(4(i)^(2) - 4(i)4(2))] 

+D pr (4^ 2 )[e feip 4(l)4(l) J RJ(2)a i + e^4(l)4(2)4(2)9 i 
+e fe ^4(l)4(l)K(2) + e Hr ^4(l)4(2)4(2)] 

+D kr {A\A 2 )Du{A\A 2 )\{L k l {l)R r ] {2) - L){l)Rl{2)){L]{l)R\{2) - 4(1)4 (2))}K (A13) 



4. WZW-term 



The last term in Eq. (fT6]) . which is proportional to the baryon number density, B°, comes from the WZW term. B° 
can be written in terms of L and R as 

B° = -^[e ij % q r(L P i (1)4(1)££(1) - 4(2)4 ? (2)^(2)) 

+ J D ps (4^2)3e^ fc ( £9rs ^(l)4(2)i 1 'I.(2) - e qrp Rt(2)L^l)Ll(l))}. (AH) 



APPENDIX B: SPIN-ISOSPIN PROJECTION 



In this appendix, we show the procedure of spin-isospin projection in detail. The first step is to replace the adjoint 
matrix in the Lagrangian with Wigner D-function. For example, 

DntiAt) -> D^iSlj.), (Bl) 

D tJ (A\A 2 ) = D mi {A 1 )D m j{A 2 ) -> D\f M , (p,x)D l MM ,^ (0 2 ). (B2) 

Next, we sandwich them between the relevant nucleon wave function, Eq . (|23|) . and integrate the Eulcr angle. We 
demonstrate the procedure by taking two examples below. For later use, we display two basic formulae for D-functions, 

J1+J2 

D* lMl (n)D* aMi (n) = Yl {JihM l M 2 \JM){J l J 2 M[M' 2 \JM')Di I ^), (B3) 
J dnD J ^* M , (Cl)Dj} 2M ,(n) = 2^1*^1. J2)S(M 1 ,M 2 )8(M{,Ml l ). (B4) 
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1. Example 1 

In this subsection, we calculate the matrix element (I 3 — J3\D\j iM , (0)|i3 J 3 ) as an example. Here and hereafter, 
we denote the nucleon state with the third component of the isospin I 3 and that of the spin J 3 simply by I/3 J3). This 
matrix element is given by the following integral, 

(h-MD\ hM ,{n)\hh) = j^i- l ) 2I:i+1 J ^Ljs^y^m^-Ui^- ( B5 ) 

Here, from Eq. (|B3[) . we note that 

3 /2 

DIum'S^D'Xj^) = E ^\M 1 -h\JM)(l\M[h\JM')D J MM ,(V) 



,7=1/2 

= (l^Mi - I 3 \\m)(i\m[J 3 \\m')D 1 ^ m ,^1) 

+{l\Mi -h\\M){l\M[Js\lM')D% 2 ul {n). (B6) 
We substitute this equation into Eq. (|B5jl and integrate the Euler angle using Eq.(jB4|. The result reads as follows, 
(h - J 3 \D\ IiM ,(n)\I 3 J 3 ) = (lhh - I 3 \lM)(l±M{J 3 \~M')S(-h,M)5(-J 3 ,M') 

= (-l) 1/2 - l3 ^S(0,M 1 )[-6(-l,Mi)5(^,J 3 ) + 5(1,M[)S(-±J 3 )], (B7) 



2. Example 2 

The second example is the matrix element of two £)-functions: (73 — J 3 \D]^ iM , (Q)D l M ^ M , (Q)\I 3 J 3 ) given by 

{h- MDl hM[ {n)Dl l2M ,{m^) 
= (2^ (_1)2/3+1 / ^-^^( n )* c U^( n ) D k^(n)^ Js (n)- (B8) 

Using Eq. (|B3p two times, we obtain 

D^^D^^D^Ij^) = (llAf 1 Af 2 |00)(HAf{M^|00') J D V / 3 J3 («) 

+ (llM 1 M 2 \\M 3 )(llM[M' 2 \lMl i ){^ 

+(i\m 3 h\\M)(l\M^-M')D% 2 MI {£l)) 

+(UM 1 M 2 1 2M 3 ) (11M[M^\ 2M^)D 2 M3M , (Q)D^ s j 3 (0). (B9) 
We substitute this equation into Eq. (|B8[) and carry out the integration of the Euler angle using Eq. (|B4|) to obtain 
(h - MD^MD^ffllhJa) = (11M 1 M 2 |1M 3 )(11A/{M^|1A^) 

(lhl 3 - I 3 \±M)(1^M' 3 J 3 \^M')5{-I 3 , M)S(-J 3 ,M') 

= (-l) 1/2 - l3 ^[6(M u l)6(M 2 ,-l) - 6(M 1) -1)6(M 2 ,1)} 
[{S(M{, 1)5{M' 2 , 0) - 8{M[,Q)S{M' 2 , l)}6{-\, h) 
-{S(M[, -l)S(M^O) + S(M[,0)d(M^ J 3 )]- (BIO) 
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